Saddle stresses for generic theories with a preferred acceleration scale 
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We show how scaling arguments may be used to generate templates for the tidal stresses around 
saddles for a vast class of MONDian theories detached from their obligations as dark matter al- 
ternatives. Such theories are to be seen simply as alternative theories of gravity with a preferred 
acceleration scale, and could be tested in the solar system by extending the LISA Pathfinder (LPF) 
mission. The constraints thus obtained may then be combined, if one wishes, with requirements 
arising from astrophysical and cosmological applications, but a clear separation of the issues is 
achieved. The central technical content of this paper is the derivation of a scaling prescription 
allowing complex numerical work to be bypassed in the generation of templates. We find that LPF 
could constrain very tightly the acceleration ao and the free parameter k present in these theories. 
As an application of our technique we also produce predictions for the moon saddle (for which a 
similar scaling argument is applicable) with the result that we recommend that it should be included 
in orbit design. 



I. INTRODUCTION 

There is an ubiquitous acceleration scale in the uni- 
verse, ao ~ 10~ 10 ms~ 2 , which turns up variously 
in cosmology and astrophysics: the cosmic expansion 
rate, galactic rotation curves, etc. This observation 
has prompted the investigation of alternative theories of 
gravity endowed with a preferred acceleration. TeVeS ll| , 
and more generally relativistic MONDian theories 043; 
provide a blueprint for such constructions. MONDian 
theories were first proposed with the motivation of by- 
passing the need for dark matter @, II • However, they 
may also be considered independently from this applica- 
tion, and be seen simply as alternative theories of grav- 
ity [9( into which an acceleration scale has been embed- 
ded. In this guise they constitute prime targets for ex- 
perimental gravitational tests inside the solar system. 

As an example let us consider TeVeS (but what fol- 
lows applies generally to what in 10] was labelled "Type 
I" theories). Abstracting from aspects which do not af- 
fect the non-relativistic limit, the theory benefits from 
the leeway of a whole free-function ^. Its choice may be 
informed by minimalism and simplicity, e.g. fj, may be 
built to encode only 2, rather than 3 or more regimes. 
Putting aside details affecting the transition between the 
two regimes, we are then left with two free parameters: 
ao (the acceleration scale of the theory) and n (control- 
ling the renormalization of the gravitational constant G). 
These are fixed by astrophysical and cosmological appli- 
cations, if the theory is to act as a competitor to dark 
matter. But ao and k can also be seen as fully free pa- 
rameters in any solar system test. 

Specifically, these theories predict a rich phenomenol- 
ogy around the saddle points of the gravitational poten- 
tial. The prospect of extending the LISA Pathfinder 
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mission so that a saddle of the Sun-Earth-Moon sys- 
tem is visited brings them within experimental strik- 
ing range [Tol - fl5j . Predictions for minimal theories con- 
strained by cosmological and astrophysical applications 
were studied in [Io[ , where the general impact of a nega- 
tive result was also examined. Detaching the target the- 
ory from its duties as "dark matter" alternative requires 
the generation of a large database of templates. However, 
re-running the adaptive-mesh code presented in [ll[ for 
each fj, is simply not feasible, and we run up against a 
computational wall. 

In this paper we show how this work can be partly al- 
leviated. As long as we are interested in changing only 
k and ao, a simple scaling argument allows the genera- 
tion of the whole set of required templates from those 
obtained with fiducial values for ao and k, shortcutting 
tedious or downright impossible hard labour. The ana- 
lytical argument is laid down in Section [Hi and its appli- 
cation to LPF is given in Section IIIII In order to include 
another topical application, in Section IIVI we also show 
how the lunar saddle would fare, were LPF to include it 
in a mission extension. 



II. SCALING BEHAVIOUR AROUND SADDLES 

Scaling is an interesting tool for generating solutions 
to apparently intractable problems. For example im- 
posing a self-similar ansatz leads to striking progress in 
the study of gravitational collapse, rendering what a pri- 
ori are PDEs into simpler ODEs (e.g. [Hill). Scaling 
behaviour was observed in the MONDian tidal stresses 
around saddles, when comparing the profiles around the 
Moon saddle and the Earth-Sun saddle (see Fig. 12 
in jllfl , and its surrounding comments) . It was noted that 
the tidal stresses are very approximately the same once 
they are spatially stretched and their amplitude scaled 
to account for the different Newtonian tidal stress A. In 
what follows we rigorously explain this empirical fact and 
extend its scope, deriving the scaling laws associated with 



2 



varying ao and n. 

In TcVeS and other type I theories [l(| the non- 
relativistic dynamics results from the joint action of the 
usual Newtonian potential $ n (associated with the met- 
ric) and a "fifth force" scalar field, (j>, responsible for 
MONDian effects. The total potential acting on non- 
relativistic particles is their sum $ = $jv + <t>- The field 
<f> is ruled by a non-linear Poisson equation: 



V • O(z)V0) = KGp, 



(1) 



with z = -r- . Here k is a dimensionless constant and 

47r ao 

ao is the usual MONDian acceleration. 

For these theories, the argument presented in Sections 
II and IV of [l2j], for a specific [i, can be generalized for 
any \x. It is always possible to define a variable: 



U = — 



Ana, 



in terms of which the vacuum equations become: 

V-U = 

a(U) U 2 V A U + U A VC/ 2 = 0, 



(2) 



(3) 
(4) 



where F is a universal function. However, we have yet 
to supply Equations (J3|) and (U]) with boundary condi- 
tions. This is done by going far enough from the saddle 
so that the field 4> has entered the Newtonian regime. 
With the conventions used in [lJJ] one has /i — > 1 (the 
renormalization in G is fully absorbed in k), and so: 



(12) 



The appropriate boundary condition is then supplied 
from the Newtonian limit relation: 
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Let us first assume that we can approximate the Newto- 
nian field around the saddle as a linear function, for the 
purpose of effectuating this matching. Then: 
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-V$iv = ArN(0,< 



(14) 



where A is the Newtonian tidal stress at the saddle point, 
and N is its angular profile (see Eqs. 35-37 in [Uj]). Defin- 
ing the MONDian "bubble size" as usual: 



i.e. the free parameters ao and k drop out, and "univer- 
sal" equations are obtained. Here 



a(U) 



dint/ 2 
din /i 



(5) 



(notice that U 2 = z 2 fi 2 (z), so that /i can be written as a 
function of U alone, /i = n(U)). The MONDian physical 
force can be obtained from U using: 



(0) 



which follows directly from @. In [12] the choice was 
made: 



fi 



(1 + C/ 2 )V4 



(7) 



so that a = 4(1 + U 2 ) 7 in agreement with Eqns. (30)-(31) 
of [TJ, and 

^-v^^(i + f / 2 )V 4 _lL, (8 ) 

in agreement with Eqn. 32 in [TJj]. However, as we see, 
the argument can be adapted to any function fi. 

Equations ([3]) and (UJ) are invariant under a rigid rescal- 
ing of the spatial variables: 



U -> U 

x — > Ax 



(9) 
(10) 



where A is spatially constant. To use the technical term 
they admit homothetic solutions, i.e.: 



U = F(Ax) 



(11) 
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(15) 



we therefore have in the Newtonian regime and close 
enough to the saddle: 



U w — N(0,< 
ro 



(16) 



This boundary condition allows us to select the homo- 
thetic solution (jlip appropriate to a given saddle and 
free parameters. To match the boundary conditions one 
should set A = I/t-q, so that the solution is 



U = F| — 

So 



(17) 



The above argument is still (approximately) valid if one 
goes beyond the linear approximation, as long as this ap- 
proximation is good up to a few ro- If the parameters ao 
and k lead to a breakdown of this assumption, however, 
then scaling is lost. 

We can now read off similar scaling laws for more fa- 
miliar quantities. Using (jSJ) we see that the MONDian 
force must have the form: 



F, = ^Gfii 



(18) 



where G is another universal function. (This scaling law 
is obvious by direct inspection of the analytical solutions 
derived for the \x used in however, as we now see, it 
is more general) . By taking derivatives we then find that 
the MONDian tidal stresses must have the form 



Sij — tzAHij ( 

Vr 



(19) 
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FIG. 1: The effect on the Signal-to-Noise Ratio (SNR) result- 
ing from varying ao, assuming different impact parameters b 
and the best estimate for the noise at the time of writing (with 
n kept fixed at re — 0.03). The fiducial value used in previous 
publications is ao = 10~ 10 ms~ 2 . Generally the larger the ao 
the higher the SNR. 



where the are also universal. This explains the scal- 
ing law observed for different A, and fixed re and ao, when 
comparing the Moon and Earth saddles But it also 
allows for templates for general values of re and oo to be 
generated from those for fiducial values simply by rescal- 
ing them according to the above laws. 

III. AN APPLICATION 

The practical applications of the previous section are 
far-reaching and will be the subject of a number of fu- 
ture publications devoted to the data analysis of a saddle 
test. As a simple example we examine in this section the 
impact of ao and re on the SNR (Signal to Noise Ratio) 
forecast for a LISA Pathfinder flyby. As in [lCj, we as- 
sume the use of an optimal noise-matched filter, using for 
our noise model the best estimate at the time of writing 
(labelled "Best Case Noise" in Fig.6 of [13]). We then in- 
spect the SNR variations with ao and re for different sad- 
dle impact parameters b. After a number of studies, fol- 
lowing on from an impact parameter b ~ 10 — 50 km 
is now considered realistic. Multiple flybys are currently 
being investigated, for which b may not be as good. We 
therefore consider SNRs for b up to 1000 km. Recall that 
for the fiducial values ao = 10 -10 ms -2 and re = 0.03 (re- 
quired, or suggested, by cosmological and astrophysical 
applications) one forecasts SNRs for the Earth-Sun sad- 
dle around 40-60 for the expected b = 10 — 50 km, only 
dropping below 5 beyond b ~ 700 km (see Fig. 7 of [HI). 

The effect of changing the acceleration scale ao is plot- 
ted in Fig. [T] It results from a change in the MOND 
bubble size ro, as predicted by Eqn. 1151 Therefore the 
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FIG. 2: Effect on the SNR obtained by varying re (keeping ao 
fixed at the fiducial value) . The fiducial value used in previous 
publications is re = 0.03. At small b, changing re may increase 
or decrease the SNR (see text for explanation). At large b one 
is better off with a small n. 
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FIG. 3: Effect on SNR obtained by jointly varying ao and 
re, for a trajectory with impact parameter b — 50km. The 
fiducial values have been indicated with a star. We also plot 
(dashed white line) the contour of constant ro passing through 
the fiducial values. 



SNR is roughly constant on lines of constant b/ao- The 
slope of the iso-SNR lines is not constant and they are 
not exactly straight because the SNR algorithm is quite 
complicated and non-linear. We see that even at large b 
it is possible to turn a weak result into a strong positive 
one by increasing ao by a factor of 2. Conversely, if ao is 
halved, a SNR below ~ 5 is now a liability for b as low 
as ~ 350 km. Without external constraints fixing ao to 
better than an order of magnitude, it is therefore risky 
to give up on a b ~ 10 — 50 km. 
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FIG. 4: Constraints placed on ao and k by a negative result for 
different impact parameters b (labelling the lines and coding 
the colours). For a given b, the admissible parameter space 
would be "outside" the corresponding b line (i.e. towards the 
right lower corner). 



FIG. 5: SNRs for the moon saddle, assuming a standard 
noise model and speed v = 0.3 km s _1 , for different impact 
parameters and day of the month (0 and 1 represent the New 
Moon, 0.5 the Full Moon). We see that the moon saddle is 
less forgiving if you miss it by more than 150 km and more 
rewarding if you get close to it (with SNRs of 200 within 
reach). If the former, new moons generate higher SNRs. 



The effect of changing k is plotted in Fig. [2] and results 
from two sources: a change in bubble size according to 
ro oc 1/k 2 (cf. Eqn. I15|) and an overall factor multiplying 
the amplitude (cf. Eqn. Q1J]). The two effects counteract 
each other, so that unless b is very large, the SNR at 
first increases with k, then decreases. For the expected 
b ~ 10 — 50 it can go either way. For large b (greater 
than b ~ 500 km for the fiducial value of ao), the bubble 
size prevails and so the SNR decreases with increasing 
k. The interplay of these two effects is best illustrated 
in Fig. [3j where we plotted the effect on the SNR of 
changing simultaneously ao and n for fixed b — 50 km. 
We also plotted the line of constant ro passing through 
the fiducial values. As we see the SNR does change along 
this line, showing that the bubble size ro is not the only 
consideration. 

Supposing we get a negative result, what constraints 
can we place upon a and re? As in [l(| we may get 
a preliminary estimate by seeking the region where the 
SNR for an optimal filter drops below 1. This was plot- 
ted Fig. @] for various values of b (in this figure, b labels 
the lines and codes the colours). For a given b, the ad- 
missible parameter space is "outside" the corresponding 
b line (i.e. towards the right-bottom corner). In gen- 
eral, a negative result forces ao to be smaller and k to be 
larger than the fiducial values, the more so, the smaller 
the impact parameter b. As we see, if we were to miss 
the saddle by 1500 km or more, the fiducial values of ao 
and k would survive a negative result. For an approach 
any closer, however, a negative result would rule them 
out and squeeze the parameter space towards the right- 
bottom corner. For 6 ~ 10 km, the ao (the k) would have 
to be smaller (larger) than the fiducial values by an order 



of magnitude. 

These constraints may now be combined with other 
pressures upon the theory, such as those arising from 
limits on G renormalisation, Big Bang Nucleosynthesis, 
fifth force solar system tests, galaxy rotation curve data, 
and cosmological structure formation. However, as advo- 
cated in the introduction, by allowing complete freedom 
in ao and k in a saddle test, we have achieved a clear 
separation of the issues confronting these theories. 



IV. THE MOON SADDLE AS A LPF TARGET 

Our technique can also be applied to a very topical 
issue: whether the Moon saddle is a good alternative 
target for LPF. Practical matters may render this saddle 
more amenable to multiple fly by s, an issue that could be 
essential in dismissing a "false alarm" , should a positive 
detection be found. In the absence of a more detailed 
study of transfer orbits we evaluate SNR's for the moon 
saddle, hoping that this may motivate further work in 
orbit design. 

Application of the algorithm in Section [TT] to the moon 
saddle is straightforward (and indeed it motivated the 
argument presented therein). As noted in [11], ro for the 
Moon saddle is smaller than the 380km found for the 
Earth-Sun saddle, and this size is more variable, depend- 
ing strongly on the phase of the Moon (it varies between 
25km and 80km; see Fig. 10 of [Hj]). However A is larger, 
too, so the tidal stresses have a larger amplitude. Nev- 
ertheless, what really matters for SNRs is the Fourier 
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transform of the signal as seen in time, with the satel- 
lite going through the bubble. The large SNRs obtained 
for the Sun-Earth saddle result from a miraculous coin- 
cidence between the sweet spot in the amplitude spectral 
density (ASD), and the size of the bubble as transformed 
into a time-signal by the typical velocities found in trans- 
fer orbits. This miracle could be spoiled by the smaller 
size of the Moon saddle. 

As it happens, orbits crossing the Moon saddle 
do so with a smaller velocity, typically smaller than 
0.5 km s" 11 . The two effects — smaller bubble, smaller 
speed — counteract each other when converting the bub- 
ble signal into a time signal. Therefore it is not surprising 
that the SNRs predicted for the Moon saddle are as high 
as those for the Earth saddle, albeit more variable in 
time, depending on the phase of the moon. 

In Fig. [5] we plotted SNRs assuming the standard noise 
model we have used throughout this paper, for a cross- 
ing of the moon saddle at v = 0.3 km s _1 , for different 
impact parameters and day of the month. On the y- 
axis and 1 represent the New Moon, and 0.5 represents 
the Full Moon. As we can see, in comparison with the 
Earth-Sun saddle, the moon saddle: 

• is less forgiving if you miss it by more than 150 km. 

• is more rewarding if you get close to it (with SNRs 
of 200 within reach). 

• in the first case, then the lunar phase is crucial, 
with the new moon producing the best results. 

In view of these results we think we should urge the orbit 
designers to include the moon saddle in their considera- 
tions. 



V. CONCLUSIONS 

To conclude, we have presented a simple argument al- 
lowing the inference of a large database of templates for 

1 Steve Kemble, private communication 



the tidal stresses that would be felt by LISA Pathfinder, 
should a saddle flyby be incorporated into the mission. 
The argument allows for the variation of the acceleration 
scale ao and k. Should the functional form of the free 
function /i be changed, the SNRs obtained would change, 
but only as predicted in |10j : they wouldn't change much 
for MONDian two-regime functions, unless b is much 
larger than ro. We may detach these theories altogether 
from their "alternative to dark matter duties" . Then we 
may consider two-regime functions with /i — > 1 at large 
arguments, but /j, oc z n , with n ^ 1, when z is small. The 
scaling argument presented here would still be applicable 
in this context, but the fiducial templates would have to 
be obtained by re-running a numerical code for each n. 
In a paper in preparation we show how this may be by- 
passed too, albeit with a much more complex analytical 
argument. 

The practical applications of our technique are far- 
reaching and will be the support of a number of future 
publications concerned with the data analysis of a saddle 
test. In this paper we merely showed how SNRs change 
by changing the parameters of the theory. This gives an 
indication of how sensitive to them the experiment is, 
and therefore how much it will constrain them. More 
importantly, as an application we applied our scaling al- 
gorithm to the prediction of results for the Moon saddle. 
The results were very encouraging and lead us to urge 
the orbit designers to include it in their considerations. 
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